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Densities of states weighted with the diagonal matrix 
elements of two operators A and B, i.e., p^^'^\E) — 

{n\A\n){n\B\n)S{E — En) cannot, in general, be written 
as a trace formula, and therefore no simple extension of semi- 
classical trace formulas is known for this case. However, from 
the high resolution analysis of quantum spectra in the semi- 
classical regime we find strong evidence that weighting the 
(5-functions in the quantum mechanical density of states with 
the product of diagonal matrix elements, (n|yl|n) (n|_B|n), is 
equivalent to weighting the periodic orbit contributions in the 
semiclassical periodic orbit sum with the product of the peri- 
odic orbit means, {A)p{B)p, of the classical observables A and 
B. Results are presented for the hydrogen atom in a magnetic 
field for both the chaotic and near-integrable regime, and for 
the circle billiard. 

PACS numbers: 05.45. -Kb, 03.65.Sq 



I. INTRODUCTION 

Semiclassical trace formulas for both chaotic |l|,|^ and 
regular Q systems relate quantum spectra and classical 
periodic orbits. These formulas have proven to be useful 
in the analysis of level statistics [Q and long-range cor- 
relations 1^ in the quantum spectra, and it has even be- 
come possible to compute individual eigenenergies from 
these expressions P-pl3[|. Gutzwiller's trace formula [^jj^] 
and the Berry- Tabor formula ^ are semiclassical approx- 
imations to the density of states but do not provide in- 
formation about experimentally measurable observables, 
i.e., matrix elements of Hermitian operators. The trace 
formulas have been extended to the calculation of diago- 
nal matrix elements of smooth operators in Refs. |pd] , |l2| . 
The extended trace formulas relate the diagonal matrix 
elements of operators to the periodic orbit means of the 
corresponding classical observables. 

However, these formulas cannot be applied directly for 
the semiclassical calculation of products of diagonal ma- 
trix element where the weighted density of states cannot, 
in general, be written as a trace formula. Products of 
diagonal matrix elements are important in several inter- 
esting applications of semiclassical theories, e.g., for the 
semiclassical theory of matrix element fluctuations [ p^ , 
with the variance of an operator A in an eigenstate |n) 



given by Var„A = {n\A^ 



(n|A|n)^. A semiclassical 



periodic orbit formula for products of diagonal matrix el- 



ements is also of crucial importance for the semiclassical 
quantization technique developed in Ref. [p^ , where the 
classical information of a set of observables is used to sig- 
nificantly improve the convergence properties of periodic 
orbit quantization. 

In this paper we investigate non-trace type formulas for 
the density of states weighted with the diagonal matrix 
elements of two operators A and B, i.e., p^^'^^{E) = 
^^ninlAln) {n\B\n)S{E — £"„). From the high resolution 
analysis of quantum spectra in the semiclassical regime 
we find strong evidence that weighting the ^-functions 
in the quantum mechanical density of states with the 
product of diagonal matrix elements, {n\A\n){n\B\n), is 
equivalent to weighting the periodic orbit contributions 
in the semiclassical periodic orbit sum with the product 
of the periodic orbit means, {A)p{B)p, of the classical 
observables A and B. 

The outline of the paper is as follows. In Sec. |l| we 
first briefly review Gutzwiller's trace formula for chaotic 
systems and the Berry- Tabor formula for integrable sys- 
tems, and discuss the extension of both equations to the 
calculation of diagonal matrix elements. We then apply 
the theories to systems with scaling properties, and intro- 
duce the high resolution analysis (harmonic inversion) of 
quantum spectra as a powerful tool to numerically verify 
the validity of the semiclassical expressions. In Sec. Ill 
we present our results on the semiclassical non-trace type 
formulas. Strong numerical evidence for the validity of 
the non-trace type equations is provided by the harmonic 
inversion of spectra of two different systems, viz. the hy- 
drogen atom in a magnetic field and the circle billiard. 
Sec. IV concludes with remarks on useful and important 
applications, and an outlook on possible generalizations 
of the non-trace type formulas. 



II. SEMICLASSICAL TRACE FORMULAS 

A. Matrix element extension of periodic orbit theory 

The quantum mechanical density of states can be 
written as the trace of the Green's function, p{E) = 
^„ S{E — En) ~ — (I/tt) Imtr G^. Replacing the quan- 
tum mechanical Green's function, = {E — H + ie)~^ 
with its semiclassical analogue and calculating integrals 
and traces in stationary phase approximation Gutzwiller 
derived the fundamental equation of periodic orbit the- 
ory i.e., the density of states expressed in terms of 



1 



quantities of the periodic orbits of the classical system. 
To obtain the density of states weighted with the diago- 
nal matrix elements of an operator A we start from the 
generalized trace formula 



P 



(A) 



(E) = - ilmtr (c+i) 

1 ^ ^ {n\A\n) 
= iim Im > 



E - En + ie 
Y,{MW)5{E ~ E^) . 



(1) 



The r.h.s. of Eq. |l| is the density of states weighted with 
the diagonal matrix elements (ri|^|7i) of the operator A. 
The semiclassical approximation to Eq. |l| for a system 
with N degrees of freedom reads pl| , |l^ 



p^^\E) = p^^\E) 



(2) 



T 



where 



.=1 y^|det(M;-/)| 



the 



term p''o'^\e) = h-^ J dqdpA{q,p)6{E - iJ(q,p)) is a 
smooth function of the energy and the fluctuating part is 
given by the periodic orbit sum, with Tp the time period, 
Sp the classical action, Mp the monodromy matrix, and 
Pp the Maslov index of the primitive periodic orbit p. 
The integer r is the repetition number of the orbit. The 
weights Ap in the periodic orbit sum are the means 
of the observable A along the periodic orbit p, i.e. 



Weyl 
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p Jo 



A{qp{t),pp{t))dt 



(3) 



The derivation of Eq. ^ requires smoothness of the ob- 
servable A over regions in phase space of size h'^ [ p^ . 
A rigorous mathematical proof of the semiclassical trace 
formula (B) using a coherent states decomposition can be 
found in |[l5]. 

In Rcfs. I ll|,|l^,|l5| formulas for the semiclassical cal- 
culation of diagonal matrix elements are obtained for 
chaotic systems with isolated periodic orbits. For regular 
systems the semiclassical trace formula for the density 
of states has been derived by Berry and Tabor For 
simplicity we restrict ourselves to systems with two de- 
grees of freedom. Assuming now that the Hamiltonian is 
classically integrable, one can express it in action-angle 
variables (I, (p) with ipi, ip2 G [0, 2tt] as H{T). For a given 
torus, LJi — dH/dIi {i — 1,2) are the corresponding an- 
gular frequencies. Periodic orbits are associated with tori 
such that the rotation number a = uji/lu2 rational, i.e., 
a = M1/M2 with Ml and M2 integers. The fluctuating 
part of the Berry- Tabor formula reads 



— =7^- Re — 



f3/2 



pi(SM(B)/?i-f i)M-f ) 



"1 1/2 



(4) 



with M = (Mi,M2) specifying the periodic orbit, and 
Tm, <S'm, and ryM the time, action and Maslov index of 
the orbit, respectively. The function qe in is obtained 
by inverting the Hamiltonian, expressed in terms of the 
actions (/i,/2) of the corresponding torus, with respect 
to I2, viz. H{Ii,l2 = gE^h)) = E By analogy with 
Eq. ^ for chaotic systems the Berry- Tabor formula (Q) can 
now be generalized straightforwardly to the semiclassical 
calculation of diagonal matrix elements [p7[ , yielding 



p'^^\E)^p'-^\e) (5) 
^ Re ^ 



IM 



M 



,i[SM(-B)/;i-f rjM-f ] 



with 



A 



M 



27r ^27r 



dtpx / C^<^2^(^l,^2,</5l,</52) (6) 



(27r)2 7o "^7o 
the classical average of the observable A on the torus. 

B. Scaling systems 

In the following we will apply Eqs. ^ and |^ to systems 
with scaling properties. In scaling systems the classical 
phase space structure does not change for all values of an 
appropriate scaling parameter, w. The scaling parameter 
is usually some power of an external field strength or, for 
Hamiltonians with homogeneous potentials, the energy. 
Examples arc billiard systems or atoms in magnetic 
fields |p^-|20|. In scaling systems the shape of periodic 
orbits does not depend on the scaling parameter, u;, and 
the classical action scales as 



Sri 



(7) 



The scaling parameter plays the role of an inverse effec- 
tive Planck constant, i.e., w = K^^- For scaling systems 
the weighted densities of states, Eqs. || and ^ can be 
rewritten as a function of the scaling parameter w, i.e. 



^ 00 

--Rej:^,^. 



^i^|dct(M;-/)| 
for chaotic systems, and 



(8) 



M 



Af|/^K|i/2 



(9) 



for two-dimensional systems with regular dynamics. 
Note that the time periods Tp and Tm in Eqs. || and 
^ must be replaced with the scaled actions Sp and sm- 
Furthermore the time average of the classical observable 
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A (Eq. H) must be replaced with the average with respect 
to the scaled action, 



A^^ 



Jo 



-4(qp(s),Pp(s))ds 



(10) 



If an observable A is chosen which is invariant under the 
scaling of the system [or scales ^ with a constant 
exponent [3] the periodic orbit amplitudes and scaled ac- 
tions in Eqs. |8| and || do not depend on w [despite a pos- 
sible power law scaling of the amplitudes with which 
can be transfered to the l.h.s. of Eqs. ^ and ^. The 
attractive feature of scaling systems is that the semiclas- 
sical weighted density of states [or more generally the 
density of states multiplied by w^] is a superposition 
of sinusoidal functions of the scaling parameter w. The 
Fourier transforms of w^^ p'^^^ (w) should therefore exhibit 
sharp peaks at the positions of the scaled actions of the 
periodic orbits. When analyzing quantum spectra, we 
will make use of the scaling advantages in the following. 



C. Precision check of the semiclassical trace formulas 

We now wish to apply the semiclassical trace formu- 
las, Eqs. ^ and |^, to a physical system with chaotic and 
regular dynamics, respectively, and to check numerically 
the validity of the semiclassical equations. The numer- 
ical check is not motivated by doubts on the validity 
of these expressions, which have been mathematically 
proven, rather we want to introduce a powerful numerical 
technique for the high precision check of equations of this 
kind. We will demonstrate the accuracy of the method 
on the well established semiclassical trace formulas here 
and then apply the same technique to numerically verify 
our conjecture on semiclassical non-trace type formulas 
in Sec. 0. 



The semiclassical trace formulas can be tested, in prin- 
ciple, by the Fourier transform analysis of quantum spec- 
tra. The Fourier transformed spectra should exhibit 
peaks at the periods (scaled actions) of periodic orbits 
with amplitudes given by the semiclassical expressions. 
However, the transformation of spectra with finite length 
yields limited resolution only, due to the uncertainty prin- 
ciple of the Fourier transform, which implies a fundamen- 
tal restriction to high precision checks of the semiclassical 
trace formulas. We therefore adopt the method of Ref. 
[ pl[ where we introduced harmonic inversion as a high 
resolution method for the analysis of quantum spectra. 
We briefly review the basic ideas of the harmonic inver- 
sion technique and refer the reader to Ref. [^l| for more 
details. 

According to Eqs. ^ and ^ the semiclassical weighted 
density of states can be written as the sum of a smooth 
background p^"^ (w) and oscillatory modulations induced 
by the periodic orbits, 



p^^Hw)^p^^\w)+ReY,4^'>e^ 



(11) 



(A) 

The amplitudes dp and scaled actions Sp of the peri- 
odic orbits are obtained from classical calculations and 

(A) 

are in general complex quantities. The amplitudes dp 
contain the phase information determined by the Maslov 
indices of orbits and the classical means of the observable 
A given by Eqs. ^ and |l^ for regular and chaotic systems, 
respectively. Instead of using the standard Fourier anal- 
ysis to extract the amplitudes and actions, we adjust a 
finite range of the quantum spectrum by the semiclassi- 
cal expression (0) with unknown and in general complex 
parameters dp^^^nd Sp. The problem of fitting a "signal" 
piA) ^y;^ to the functional form (|l^) is known as harmonic 
inversion. As a numerical technique for the harmonic in- 
version of a signal, i.e. a quantum spectrum, we apply 
the method of filter-diagonalization |2^,^ which allows 
extracting the spectral quantities in any given interval 
of interest. Operationally, one proceeds by setting up a 
small generalized eigenvalue problem. The actions Sp in 

(A) 

the chosen spectral domain and amplitudes 

dr 

are ob- 
tained from the resulting eigenvalues and eigenvectors. 
Thus, the recurrence spectrum is effectively discretized, 
the number of terms being the number of eigenvalues in 
the spectral domain. This method is a variational one (as 
opposed to the Fourier transform) and therefore practi- 
cally has an infinite resolution once the amount of infor- 
mation contained in the signal /o'-'^'(w) is greater than 
the total number of unknowns dp"^' and Sp. 

As a physical system for the high precision analysis 
of quantum spectra and the comparison with the semi- 
classical trace formulas we choose the hydrogen atom 
in a magnetic field [nB-EQ|. This is a scaling system. 



with 



w = 7 



^/■^ = h^ff the scaling parameter and 
7 = -6/(2.35 X 10^ T) the magnetic field strength in 
atomic units. Introducing scaled coordinates ^■^^^r and 
momenta 7~^/'^p and choosing the projection of the an- 
gular momentum on the magnetic field axis = one 
arrives at the scaled Hamiltonian 



(12) 



with E = E^^"^/'^ the scaled energy. The classical dy- 
namics is near-integrable at low energies, E < —0.5, and 
undergoes a transition from regularity to chaos in the 
energy region —0.5 < E < —0.13. At energies above 
E = —0.13 a Poincare surface of section analysis of the 
classical dynamics does not exhibit any regular structures 
larger than of microscopic size [|l^ . We compare spectra 
at constant scaled energy E = —0.1 with the results of 
the semiclassical trace formula (^) for chaotic systems, 
and spectra in the near-integrable regime at i? = —0.5 
with the extended Berry- Tabor formula (^. We choose 
two different operators. The first. 



A = 



1 



2rp2 



(13) 



has already served to study the distribution of transition 
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matrix elements in classically chaotic and mixed quan- 
tum systems P4LE5|. The second operator is 



B = rp^ 



(14) 



Eigenvalues of the scaling parameter w are obtained by 
solving Schrodinger's equation (in semiparabolic coordi- 
nates /i — \Jr \- z and v = — z) 

1 



with the radial operators pj^ and defined as 



(15) 



1 d 



d 



.2 1 d 
Pu = — ^ I 



Eq. can be written in matrix form by using an ap- 
propriate basis set. The resulting generalized eigenvalue 
problem is solved numerically. It has to be noted that 
the eigenvectors obtained, IV'n), are orthonormal with 
respect to the scaled momentum operator, i.e. 



{'^Pyn\pi + Pl\lpn) = {m\n) = (5„ 



(16) 



with modified eigenvectors \n) defined by \n) = (p^ + 
p2)i/2|^^) The diagonal matrix elements of an operator 
A are therefore obtained as 



{n\A\n) = {i,^\A{pl+pl)\'il,^ 



(17) 



We are now prepared to compare the quantum spectra of 
the hydrogen atom in a magnetic field with the semiclas- 
sical approximations in the chaotic and regular regime of 
the classical phase space. 



1. Chaotic regime 

We have calculated 3181 eigenvalues Wn < 80 of the 
scaling parameter and the diagonal matrix elements of 
the two operators A — l/(2rp2) and B = rp^ for the hy- 
drogen atom in a magnetic field at constant scaled energy 
E = —0.1. The distributions of the matrix elements are 
presented in Fig. |l|. The matrix elements are distributed 
randomly around the mean values without showing any 
regular pattern, as is typical of systems with chaotic dy- 
namics. The quantum mechanical weighted density of 
states 



P 



(A), 



w) — ''^^{n\A\n)S{w — Wn) 



(18) 



can now be analyzed with the harmonic inversion tech- 
nique to obtain the scaled actions Sp and the amplitudes 

d^^^ (see Eq. ^l]) of the classical periodic orbits. As can 
be seen from Eq. the periodic orbit amplitudes 
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FIG. 1. Values of the diagonal matrix elements (n|A|n) 
for the hydrogen atom in a magnetic field at scaled energy 
E — —0.1 in the chaotic region of phase space as functions 
of the scaling parameter w; = 7 
strength): (a) operator A — l/(2rp'^ 



(7 ~ magnetic field 
(b) B = rp^ 



A.pdp 



(19) 



are given as the product of the amplitudes, dp of 
Gutzwiller's original trace formula, and the classical pe- 
riodic orbit means Ap in Eq. |l^. For the graphical pre- 
sentation of the results it is therefore convenient to di- 
vide the quantum amplitudes dp^^ obtained by the har- 
monic inversion of the spectra by the amplitudes, dp of 
Gutzwiller's trace formula. The periodic orbit quantities 
Ap obtained in this way from the quantum spectra at 
scaled energy E = —0.1 are presented in Fig. || for three 
different operators, viz. (a) the identity /, and the oper- 
ators (b) A — l/(2rp2) and (c) B = rp^ . The solid lines 
and crosses in Fig. |^ mark the periodic orbit means ob- 
tained by the harmonic inversion of the quantum spectra. 
For comparison the dashed lines and squares present the 
periodic orbit means of the observable obtained classi- 
cally via Eq. 10, For the identity the classical periodic or- 
bit averages (squares in Fig. ||a) are exactly equal to one. 
This is in excellent agreement with the harmonic inver- 
sion analysis of the quantum mechanical density of states 
(crosses in Fig. ^), despite the two weakly separated pe- 
riodic orbit contributions around s/2tt « 1.1. For the 
two operators A = l/(2rp2) in Fig. ||b and B = rp^ in 
Fig. Ep the agreement between the periodic orbit means 
obtained by harmonic inversion of the quantum spectra 
and classically by Eq. is of similar high accuracy as 
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FIG. 2. Periodic orbit means of observables (a) the iden- 
tity, (b) A = l/(2rp ), and (c) B — rp for the hydrogen 
atom in a magnetic field at scaled energy E — —0.1 as func- 
tions of the scaled action. Solid lines and crosses: Results of 
the harmonic inversion of quantum spectra. Dashed lines and 
squares: Periodic orbit means obtained by classical calcula- 
tions. The agreement between the quantum and the classical 
calculations seems to be excellent, except for the nearly de- 
generate recurrences at s/2n ~ 1.1. 

for the identity in Fig. ||a. The results presented in Fig. |^ 
demonstrate that harmonic inversion of quantum spectra 
pl[ is indeed a powerful tool for the high precision check 
of semiclassical theories. Fig. || provides an excellent nu- 
merical verification, by way of example of the hydrogen 
atom in a magnetic field and the chosen set of operators, 
of the validity of the semiclassical trace formula (|^) for 
chaotic systems. 



2. Regular regime 

In the same way as described above we have checked 
the validity of the extended Berry- Tabor formula (^) for 
integrable systems. As a physical system we again choose 
the hydrogen atom in a magnetic field, but at low scaled 
energy E = —0.5, where the classical phase space is reg- 
ular. We have calculated 5640 eigenvalues w„ < 160 of 
the scaling parameter and the diagonal matrix elements 



in the same way as explained above. The results obtained 
for the regular system at scaled energy E = —0.5 resem- 
ble those of Fig. || for the chaotic system. The difference 
is that the averages of the observables for the resonant 
tori have been extracted from the quantum spectra by 
application of the generalized Berry- Tabor formula (||). 
The quantum results perfectly agree with the classical 
averages which illustrates the validity of the generalized 
Berry- Tabor formula. 



III. NON-TRACE TYPE FORMULAS 

The generalized semiclassical trace formulas (||) and 
discussed in Sec. || allow the semiclassical calcu- 
lation of the diagonal matrix elements of smooth op- 
erators. However, it would be desirable to know even 
more generalized expressions for the calculation of prod- 
ucts of matrix elements. As mentioned in the introduc- 
tion, such formulas arc important, e.g., in the semiclas- 
sical theory of matrix element fluctuations ||l^ or for the 
construction of cross-correlated periodic orbit sums . 
To study matrix element fluctuations of an operator A 
the density of states can be weighted with the variances 
Vai„A = - {n\A\ny, i.e. 

^(Var ^ J2{n\A'\n)SiE - E^) 

n 

-J2{n\MnfS{E~E^). (20) 

n 

The first term in (GO) can be written as a semiclassical 
trace formula (Eqsrj2 and |^ for chaotic and regular sys- 
tems, respectively) with the observable A replaced with 
its square, A^. However, because of the squares of the 
matrix elements, the second term in ( pO| ) in general can- 
not be expressed in a straightforward fashion with the 
help of the Green's operator as a trace formula. The 
trivial exception is when the operator A commutes with 
the Hamiltonian, which means that A is a constant of 
motion and thus its variance vanishes. Thus the deriva- 
tion of a semiclassical approximation to the second term 
in ( pO| ) constitutes a nontrivial problem. 

One solution can be obtained by application of periodic 
orbit sum rules [Q. Using smooth approximations of the 
(5-functions, e.g. Gaussians of width e, 



S.iE) = 



1 



and the relation 



27re 



(21) 



(22) 



of the two operators A = 

weighted density of states ( |18| ) for the identity, and the 
operators A = l/(2rp^) and 13 = rp^ have been analyzed 



2rp^) and B = rp\ The the second term in iM) can formally be written as the 



square of the densityof states weighted with the diagonal 
matrix elements 10,131, viz. 
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Y,{n\A\n)H^^^,{E-Er.) 

71 

2y^eJ2i^\Mnf^eiE - En 



(23) 



The width e in (|23| ) must be chosen sufficiently small so 
that the smoothed (5-functions do not overlap. However, 
it should be noted that this condition cannot be fulfilled 
for systems with degenerate states. On the r.h.s. of Eq. 
p3| the weighted density of states can now be written as a 
trace formula and replaced with its semiclassical approx- 
imations (H) and (H) for chaotic and regular systems, re- 
spectively. Evaluation of the square of the periodic orbit 
sum on the r.h.s. of ( p3| ) then finally yields a double sum 
over all periodic orbits of the classical system. Although 
this result is formally correct, it is very inconvenient for 
practical applications for the following reasons. Firstly, 
the number of periodic orbits proliferates exponentially 
in chaotic systems and the handling of the single periodic 
orbit sum is already nontrivial. The practical evaluation 
of the double sum would be even more cumbersome. Sec- 
ondly, the width e in ( p3|) is a free parameter. Although 
the results should not depend on the width if e is chosen 
sufhciently small, the appropriate choice may render nu- 
merical calculations extremely expensive. Thirdly, and 
most importantly, the r.h.s. of Eq. ^ does not coincide 
with the "simple" trace formulas in those special cases, 
where the operator A commutes with the Hamiltonian. 
Even for the simplest operator, the identity A = /, we 
end up with the nontrivial periodic orbit sum rule of Ref. 
Q instead of Gutzwiller's trace formula for the density of 
states. Especially the third point indicates that the pro- 
cedure described above might not be the simplest way 
to construct a semiclassical approximation to non-trace 
type formulas such as Eq. |2^. It is the main objective of 
this section to present a semiclassical approximation to 
non-trace type weighted densities of states. Our semiclas- 
sical expressions agree with the well established "simple" 
semiclassical trace formulas when the weighted density 
of states can be written, for at least one of the operators 
commuting with the Hamiltonian, as a quantum mechan- 
ical trace formula. 

Starting from a more general equation than (|2C 
study the density of states 



P 



{A,B) 



{E) = ^{n\A\n){n\I3\n)5{E - E^ 



(24) 



weighted with the product of the diagonal matrix el- 
ements of two smooth operators A and B. Eq. |2j 
is the starting point to construct a quantum mechan- 
ical cross-correlation function from a set of operators 
m Ref. g. The variance of matrix elements (Eq. pO| ) 
is obtained by setting B = A. The weighted density 
of states (24) can only be written as a trace formula. 



P^^'B\E) = (-l/7r)Im iT{AG%B} if either i or B com- 
mutes with the Hamiltonian. As discussed in Sec. || (see 
Eq. |l^ the semiclassical expressions for the weighted 
densities of states differ from Gutzwiller's trace formula 
and the Berry- Tabor formula in the following way. The 
periodic orbit amplitudes are multiplied with the classi- 
cal periodic orbit (or torus) averages of the observable 
A. We now assume that this ansatz is still valid for the 
non-trace type weighted density of states (p^, i.e., its 
semiclassical analogue has the same functional form as 
Gutzwiller's periodic orbit sum but with periodic orbit 
amplitudes dp multiplied with the classical averages Ap 



and Bp of both observables A and B, 



(25) 



with Ap and Bp given by Eqs. ^ and g for chaotic and 
regular systems, respectively. As can easily be seen, this 
ansatz has the property that the trace formulas (||) and 
(^ are recovered if one of the operators is chosen to be 
the identity or one of the operators commutes with the 
Hamiltonian. However, the general validity of this ansatz 
is not at all obvious, and will be checked numerically by 
the high resolution analysis of quantum spectra in the 
following. With the ansatz (|2^) for the periodic orbit 
amplitudes the semiclassical analogue to the non-trace 
type formula (M) reads 



I oo 

^ReJ^ApBpY: 



^1 y^|det(M;-/)| 



(26) 

JlSp{E)/h-^lip]r 



for systems with underlying chaotic classical dynamics, 
and 



P^^^^He) = p[,^'^\E) 



(27) 



^ Re V A^B^ ^ ^.[sME)/n-^v^-^] 



M 



for integrable systems. Eqs. ^ and [2^ are the central 
propositions of this paper, and generalize the semiclas- 
sical trace formulas (||) and (0) to the non-trace type 
weighted density of states (P^. The nontrivial state- 
ment is that weighting the quantum mechanical density 
of states with the product of diagonal matrix elements 
of smooth operators is equivalent, on the semiclassical 
level, to weighting the periodic orbit contributions in the 
periodic orbit sum with the product of the averages of 
the corresponding classical observables. 

In analogy to the discussion of scaling properties in 
Sec. H, Eqs. |6| and ^ can be reformulated for scahng 
systems, viz. 



^ OQ 



(28) 



r^i ^\dei{M;-I)\ 
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for chaotic systems, and 



(29) 



Rc ^ Am^m- 



M 



,'i[s]Vi'«-f »)lV[- 



''|l/2 



for two-dimensional systems with regular dynamics. For 
scaling systems the classical periodic orbit averages Ap 
and Bp in ( p^ ) must be calculated with respect to the 
classical action instead of time as defined in Eq. 

In the following we will provide convincing numerical 
evidence for the validity of the semiclassical non-trace 
type formulas by the high precision analysis (harmonic 
inversion) of quantum spectra of two different systems, 
viz. the hydrogen atom in a magnetic field and the circle 
billiard. A rigorous mathematical proof of the expres- 
sions given above is still lacking and constitutes a chal- 
lenge for the further development of semiclassical theo- 
ries. 



A. Hydrogen atom in a magnetic field 

To demonstrate the validity of the semiclassical non- 
trace type formulas, Eqs. and we use the same sys- 
tem and set of operators as in Sec. viz. the hydrogen 
atom in a magnetic field at scaled energies E — —0.1 and 
E = —0.5 in the chaotic and near-integrable regime, re- 
spectively, and the operators A = l/(2rp^) and B — rp^. 
With the quantum mechanical eigenvalues and diagonal 
matrix elements at hand, we construct the weighted den- 
sities of states (see Eq. ||) (a) p'^'^'^^w), (b) p<-^^^^w), 
and (c) p''^'^\'w). These spectra are analyzed with the 
harmonic inversion technique as described in Sec. The 
analysis provides the scaled action Sp of the periodic 

orbits and the periodic orbit amplitudes d^'^^ [Sj^'^^ 
and d^'^''). The amplitudes of the weighted densities of 
states are divided by the amplitudes of the unweighted 
densities of states to obtain, according to Eq. pa, the 
products of the periodic orbit means Ap [Bp and ApBp). 
These values are presented as solid lines and crosses in 
Fig. ^ for the spectra in the chaotic regime at scaled 
energy E = —0.1 and in Fig. ^ for the spectra at scaled 
energy E — —0.5 in the near-integrable regime. For com- 
parison, the squares mark the products of the periodic 
orbit means obtained from the classical calculations. As 
in Fig. H for the high precision check of the semiclassical 
trace formula (H), the agreement between the quantum 
and classical calculations is found to be very good, which 
strongly supports the validity of the semiclassical non- 
trace type expressions. Note that the somewhat larger 
deviations between the crosses and squares for the nearly 
degenerate recurrencies at Sp/2'K w 1.1 in Fig. ||have also 
been observed in Fig. ^ for the semiclassical trace formu- 
las, i.e., the deviation does not indicate any failure of the 
non-trace type formula (|2 
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FIG. 3. Products of periodic orbit means of the two ob- 
servables A = l/(2rp'^) and B = rp^ for the hydrogen atom 
in a magnetic field at scaled energy E = —0.1 as functions 
of the scaled action. Solid lines and crosses: Results of the 
harmonic inversion of the non-trace type weighted densities of 
states. Dashed lines and squares: Results obtained by clas- 
sical calculations. As in Fig. ^, the agreement between the 
quantum and the classical calculations seems to be excellent, 
except for the nearly degenerate recurrences at s/2ti ~ 1.1. 



B. Circle billiard 

We now investigate the validity of the semiclassical 
non-trace type formula (^9|) on a second system, viz. 
the integrable circle billiard. This system also serves 
as a model example in Ref. to construct a semi- 
classical cross-correlated periodic orbit sum for a given 
set of smooth observables, and to calculate semiclassical 
spectra and diagonal matrix elements by harmonic inver- 
sion of the cross-correlation function. As is well known, 
Schrodinger's equation for the circle billiard with radius 
R can be separated in polar coordinates (r, Lp) , and the 
wave functions can be expressed in terms of Bessel func- 
tions. 



V'nm(f, ^fnmJ\m\iknmr)e 



(30) 



with the J\fnm being normalization constants, m the 
angular momentum quantum number, and knm = 
\/2M Enm/fi' the quantized wave numbers obtained as 
the nth zero of Bessel functions, J|„j| (A:„mi?) = 0. In 
the following we choose radius R = 1. We calculated 
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FIG. 4. Same as Fig. 5 but at scaled energy E — —0.5. 

31208 eigenvalues knm < 500, and the diagonal matrix 
elements of the operators r and r^. The quantum spec- 
tra of (a) the unweighted density of states p{k), (the 
wave number k is the scaling parameter, w — k for bil- 
liard systems 113]) and the density of states weighted 



with the matrix element expressions (b) 



(c) i^nmlrlipnm) , and (d) the variance (Var r) 
(V'nmk^lV'nm) " (V'nm |r|i/'nm) ^ havc been analyzed with 
the harmonic inversion method. The amplitudes ob- 
tained, divided by the amplitudes of the Berry- Tabor 
formula, are presented as solid lines and crosses in Fig. 
^, and the corresponding classical averages are drawn 
as squares for comparison. As can be seen, the agree- 
ment is perfect, not only for the identity and the pe- 
riodic orbit means of the observable in Fig. ^ and 
verifying the Berry- Tabor formula and its exten- 
sion (^, but also for the squares of the periodic orbit 
means of r and the variance of this observable in Fig. 
^ and where the agreement demonstrates the va- 
lidity of the non-trace type equation (^) for the cir- 
cle billiard with A — B — r. The squares in Fig. ||d 
mark the classical variances of the observable r on the 
various resonant tori. Our conjecture therefore provides 
a basic formula for semiclassical matrix element fluctu- 
ations, since it directly relates the quantum variances 
Var„^ = — of a smooth operator A 

to the classical variances Varp^l = {A'^)p — {A)^ of the 
observable A taken along the periodic orbits or resonant 
tori. 

The perfect agreement between the quantum and clas- 
sical results for the circle billiard in Fig. || compared to 
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FIG. 5. 
the identity. 

Varpr = (r^)p — (r)p for the circle billiard with radius R = 1 
as functions of the scaled action. Solid lines and crosses: Re- 
sults of the harmonic inversion of quantum spectra. Dashed 
lines and squares: Periodic orbit means obtained by classical 
calculations. 

the very good but not absolutely perfect results for the 
hydrogen atom in a magnetic field in Figs. ^ and ^ may be 
explained by the different number of quantum states used 
in the harmonic inversion analysis. For the circle billiard 
we have calculated more than 30000 states, which is by 
about a factor of 10 (5.5) times more quantum states 
than for the hydrogen atom in a magnetic field at scaled 
energy E = -OA {E = -0.5). 



IV. CONCLUSION AND OUTLOOK 

We have extended semiclassical trace formulas for the 
density of states of regular and chaotic systems, or the 
density of states weighted with the diagonal matrix ele- 
ments of smooth operators, to the more general class of 
non-trace type equations, where the density of states is 
weighted with the diagonal matrix elements of two opera- 
tors A and B, i.e., p^^'^^E) = J2n{n\Mn){n\B\n)5{E - 
En). By the high resolution analysis (harmonic inver- 
sion) of the quantum spectra of two different systems, 
viz. the hydrogen atom in a magnetic field and the circle 
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billiard, we have given numerical evidence that weighting 
the quantum mechanical density of states with the prod- 
uct of the diagonal matrix elements {n\A\n){n\B\n) is 
equivalent, on the semiclassical level, to weighting the pe- 
riodic orbit contributions in the periodic orbit sum with 
the product of the averages of the corresponding classical 
observables, {A)p{B)p, where the means are taken along 
the periodic orbits or resonant tori for chaotic and regular 
systems, respectively. However, a rigorous mathematical 
derivation of semiclassical non-trace type formulas ap- 
pears nontrivial, and would be a challenging task for the 
further development of semiclassical theories. 

There are several useful and important applications 
of semiclassical non-trace type formulas. For example, 
it enables the semiclassical approach to matrix element 
fluctuations. The variances Var„^ = (n|yl^|n) — (n|A|n)^ 
of the diagonal matrix elements of a smooth operator 
A are expressed in terms of the variances YaipA = 
{A^)p — {A)p of the classical observable A taken along 
the periodic orbits or resonant tori. Non-trace type for- 
mulas also provide the semiclassical approximation to 
cross-correlated weighted density of states, Paa' (E) = 
J^niM^aln) {n\Aa'\n)S{E — En) with a set of smooth op- 
erators Aq, a = 1, . . . , D. The additional classical infor- 
mation obtained from the set of classical observables can 
be used to significantly improve the convergence proper- 
ties of semiclassical quantization methods |14| . 

In this paper we have investigated non-trace type ex- 
pressions for products of two diagonal matrix elements. 
These products have been chosen because of the impor- 
tant applications to semiclassical matrix element fluctu- 
ations, i.e., the calculation of variances of matrix ele- 
ments and to the semiclassical quantization method in 
Ref. iQ. However, our conjecture is not restricted to 
products of two matrix elements. For example, it ap- 
pears straightforward to generalize Eqs. and to 
products of more than two matrix elements and classi- 
cal periodic orbit means. The most general case of non- 
trace type equations would be the analysis of functions 
f(Ann, Bnn, C'nn , • ■ ■) of One Or morc diagonal matrix el- 
ements, i.e., p'^^^E) = J2nf('^nn,Bnn,Cnn,---)S{E - 

En), which should be obtained semiclassically by multi- 
plying the periodic orbit amplitudes of Gutzwillcr's trace 
formula or the Berry- Tabor formula with the function 
f{{A)p, {B)p, {C)p, . . .) of the periodic orbit means of the 
corresponding classical observables. Certainly the opera- 
tors and the function / must be smooth. Clearly, further 
investigations will be necessary to verify that conjecture 
and to specify the smoothness conditions on operators 
and functions. 

In conclusion, the analysis of non-trace type equations 
will provide a valuable instrument for extending the re- 
lation between quantum mechanical matrix elements on 
the one side and the periodic orbit means of classical 
observables on the other. 
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